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ABSTRACT 

We investigate the dynamics of dust grains with various sizes in protoplanetary disk winds driven 
by magnetorotational turbulence, by simulating the time evolution of the dust grain distribution in 
the vertical direction. Small dust grains, which are well coupled to the gas, are dragged upward with 
the up-flowing gas, while large grains remain near the midplane of a disk. Intermediate-size grains 
float near the sonic point of the disk wind located at several scale heights from the midplane, where 
the grains are loosely coupled to the background gas. For the minimum mass solar nebula at 1 AU, 
dust grains with size of 25 - 45 iim float around 4 scale heights from the midplane. Considering the 
dependence on the distance from the central star, smaller-size grains remain only in an outer region of 
the disk, while larger-size grains are distributed in a broader region. We also discuss the implication 
of our result to the observation of dusty material around young stellar objects. 

Subject headings: accretion, accretion disks - dust, extinction - ISM: jets and outflows - planets and 
satellites: formation - protoplanetary disks - stars: pre-main sequence 


1. INTRODUCTION 

Protoplanetary disks are believed to be the birth place 
of planets (e.g., Safronov 1972; Hayashi et al. 1976, 1985; 
Cameron 1978), and have been extensively studied from 
both theoretical and observational sides (e.g., Armitage 
2011, and references therein). One of the major un¬ 
knowns in the evolution of protoplanetary disks is the 
dispersal mechanism at the Hnal stage. From near in¬ 
frared observations the typical dissipation time of pro¬ 
toplanetary disks is inferred to be 10® — 10^ yr (Haisch 
et al. 2001; Hernandez et al. 2008; Yasui et al. 2009; Tak- 
agi et al. 2014, 2015). The photoevaporating wind by 
UV(e.g., Hollenbach et al. 1994; Alexander et al. 2006) 
and X-ray (e.g. Ercolano et al. 2008) emissions from 
a central star has been extensively investigated. The 
disk wind driven by magneto-turbulence, which is the 
magneto-rotational instability (MRI hereafter; Velikhov 
1959; Chandrasekhar 1961; Balbus & Hawley 1991) in 
a disk, is also proposed as a cooperating process or a 
powerful alternative in dispersing the gaseous component 
(Suzuki & Inutsuka 2009). 

In order to promote the mass accretion by the trans¬ 
port of the angular momentum in a protoplanetary disk 
(Lynden-Bell & Pringle 1974), the gas has to be in a tur¬ 
bulent state at least in some portions of the disk to keep 
the effective turbulent viscosity since the microscopic vis¬ 
cosity is negligible. Among various mechanisms, MRI is 
supposed to be a promising candidate to excite such tur¬ 
bulence. Suzuki & Inutsuka (2009) showed that the MRI- 
driven turbulence inevitably induces vertical outflows by 
the Poynting flux of the magnetohydrodynamical tur¬ 
bulence (see also Miller & Stone 2000), and such disk 
winds potentially contribute to the dispersal of the gas 
component of protoplanetary disks (Suzuki et al. 2010), 
whereas this processes has been extended (Ogilvie 2012; 
Bai & Stone 2013a; Lesur et al. 2013) and the quanti- 
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tative mass flux carried by the disk wind is still under 
debate (Fromang et al. 2013). 

A spectro-astrometric survey of protoplanetary disks 
by Pontoppidan et al. (2011) showed non-Keplerian mo¬ 
tions that are consistent with wide-angle disk winds, in 
addition to the component of Keplerian rotation. This 
might be an indirect evidence of such MRI turbulence- 
driven vertical outflows. Natta et al. (2014) also reported 
the observations of slow gaseous outflows with a few 
km s“^ from protoplanetary disks around T-Tauri stars. 
Ellerbroek et al. (2014) observed time-variable optical 
fadings and infrared brightenings, which is supposed to 
reflect dust grains lifted up high above the disk midplane 
by a slow wind. Measurement of hydrogen atoms toward 
young stars show the absorption of the Lyman-a band in 
excess of that is inferred from the circumstellar extinc¬ 
tion (McJunkin et al. 2014), which may be also caused by 
dust grains stirred up to the upper layers. Near-infrared 
imaging of RY Tau suggested that the disk could have 
an optically thin and geometrically thick layer above the 
disk surface (Takami et al. 2013), which may reflect dust 
grains floating at the upper atmosphere. 

The dynamics of dust grains in the MRI turbulence- 
driven disk winds should be potentially interesting in 
the context of planetesimal formation but have not been 
studied yet. In this paper, we investigate the time evolu¬ 
tion of the vertical distribution of dust grains with vari¬ 
ous sizes in the disk wind by simple one-dimensional (ID) 
simulations. 

2. MODEL 

We compute the time evolution of the vertical profile 
of dust density (Section 2.1) under a fixed background 
gas component (Section 2.2) in a (ID) vertical box from 
z = —lOiLo to z = lOiLo, where Hq = is the 

pressure scale height derived from the sound speed, Cg, 
and the Keplerian rotation frequency, Hk. In order to 
avoid unphysical effects at the boundaries, we set up an 
extra ±5i7o region outside the top and bottom boundary. 
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and at z = zLIBHq the outflow boundary condition is 
prescribed. As the initial condition, all the dust grains 
are sedimented at the midplane, z = 0, and the total 
mass of the dust in the computational domain is 1/100 
of the total mass of the gas. We use the normalization 
of ilk = 1 and Hq = 1 ) which gives Cg = 1/2. 


2.1. Dust Component 

Dust grains are treated as “pressure-less fluid” in a 
fixed background gas distribution that is modeled in Sec¬ 
tion 2.2, namely drag force is exerted on the dust fluid 
from the background gas component, but the back reac¬ 
tion from the dust to the gas is not taken into account. In 
this work, we consider dust grains of which size is smaller 
than a mean free path of the gas. In this case, the gas 
drag force is in the Epstein’s regime (Adachi et al. 1976; 
Sekiya 1983; Paardekooper 2007). For example, in the 
minimum mass solar nebula (MMSN hereafter; Hayashi 
1981), the mean free path is > 1 cm for distance > 1 AU 
from the central star (Nakagawa et al. 1986); the dust 
grains we consider here is smaller than this value. Then, 
the drag force is written as (Schaaf 1963) 

Fd = —n pga^\l\v\l\v 


where pg is the gas density, a is the size of a dust 

grain, = (8/7r)^^^Cs is the mean thermal velocity 
(Takeuchi & Lin 2002), Av is the relative velocity of 
dust to gas and A4 is Mach number of the relative flow, 
A4 = |Au|/cs. Equation (1) has the following asymptotic 
behavior (Adachi et al. 1976; Paardekooper 2007), 

F - I if |Au| < Cs 

^ \ —7rpga^|Au|Au if |Au| ^ Cg. ^ 


In the small Av limit, the drag force is from stochas¬ 
tic collision by randomly moving gas particles, and thus, 
the collision rate is proportional to the thermal veloc¬ 
ity, which gives the dependence of Fd on vth- On the 
other hand, in the large An limit, the dominant factor 
that determines the collision rate is the relative veloc¬ 
ity between dust and gas. Therefore, the dependence on 
nth in the drag force for the small An limit is replaced 
by the dependence on An for large An. In this work, 
instead of using Equation (I), we adopt a simplified for¬ 
mula that interpolates between these limits (Kwok 1975; 
Paardekooper 2007), 

47 r n 

-F’d = — -pga^vthfdlAvl, (3) 

where _ 

The stopping time for the gas drag is defined as tgtop = 
mAv/ |Fd|, where m is the mass of a dust grain, tgtop 
is a time scale for the decay of the relative velocity be¬ 
tween the dust and the gas. By normalizing with 11x7 
define the nondimensional stopping time, Tg = fgtopf^k- 
In this paper we assume that the dust grains have spher¬ 
ical shape with homogeneous internal density, pm- Since 
the mass of a dust grain is m = (4/3)7ra^Pm7 then the 


nondimensional stopping time is rewritten as 


Ts 


P^a 

PgVth 




(5) 


For later discussion, we evaluate Tg at the midplane 
of the MMSN. The relative velocity of dust to gas is 
much smaller than the sound speed around the midplane. 
Therefore, the nondimensional stopping time, Tg^mid? at 
the midplane is 



where we adopt pm = 2 g cm“^ in this work (Youdin 
& Lithwick 2007), and subscripts “mid” and “lAU” 
denote that the values are evaluated at the midplane 
and at r = 1 AU, respectively. Here, we use, in 
addition to the relation for the Keplerian rotation, 
Dk = Hk,iau(?'/ 1AU)“^/^, the scalings for the den¬ 
sity at the midplane and the sound velocity of the 
MMSN (Hayashi 1981; Sano et al. 2000): Pg.mid = 

Pg.mid.iAU (»"/lAU)““^^ with Pg,mid.iAU = 1-4 X 10“®g 
cm“^ and Cg = Cg^AU (r/lAU) with Cg^AU = 0.99 
km s“^. The nondimensional stopping time, Tg, mea¬ 
sures whether dust grains are tightly (ts < 1) or weakly 
(Tg > 1) coupled to the background gaseous flow. 

The velocity of dust fluid is subject to the following 
momentum equation (Takeuchi & Lin 2002): 


dvd 

dt 




J 


tstop Pdtgtop 


( 7 ) 


where Vd and Vg are the vertical velocity of dust and gas. 
The first term is the vertical component of the gravita¬ 
tional force by the central star, the second term denotes 
the gas drag force on the dust from the mean gas flow, 
and the third term represents turbulent diffusion. J is 
the mass flux due to turbulent diffusion and modeled as 


Pg^ d (Pd\ 
Sc dz \pgj 


( 8 ) 


(Takeuchi & Lin 2002), where n is effective kinetic vis¬ 
cosity of gas and expressed as n = ai7QHk/2 with a 
parameter (Shakura & Sunyaev 1973). Here, we adopt 
a constant a = 0.01 in the entire calculation region for 
simplicity. The value of a controls the timescale of the 
upward diffusion of grains that are located initially at 
2 : = 0 to the wind region, 2 ; > SFq. However, the final 
profile of pd in the quasi-steady state does not depend 
on a provided that a is finite, although smaller a re¬ 
quires longer time to achieve the quasi-steady state. The 
Schmidt number Sc is a ratio of gas and dust diffusion 
coefficients, DgjD^, which represents the strength of the 
coupling between the gas and the dust. For small grains. 
Sc approaches unity so that the grains have the same 
diffusivity as the gas, while it becomes infinite for large 
grains. For intermediate grains, it is S'c « l + (Youdin 
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& Lithwick 2007). In this paper, we use Sc = 1, because 
smaller grains with Ts < 1 are tightly coupled to the 
gas. Therefore, the steady state condition is achieved for 
the dust, and the dust velocity reaches terminal veloc¬ 
ity. Thus, applying the time-steady condition (^ = 0) 
to Equation (7), we have 

Vd= Vg- V,kTsZ + —. (9) 

Pd 

In addition, by substituting Equation (9) into the conti¬ 
nuity equation of dust, ^ -I- V • (pd'i’d) = 0, we obtain 

^ ^ (PdFd + J) = 0, (10) 

where we separate the dust velocity Vd into the mean 
flow, 

Vd ~ W (H) 

and the diffusion component, J/pd- 

By solving Equation (10), we follow the time evolution 
of the vertical distribution of the dust fluid. 

2.2. Gas Component 

We adopt the result of the three-dimensional (3D) 
MHD simulation in a local shearing box by Suzuki & In- 
utsuka (2009) for the fixed background gas component. 
In order to incorporate it into our vertical ID calcula¬ 
tion, the time- and horizontal plane-average is taken for 

Vg. 

Since our computation region, \z\ < IOHq, is larger 
than the region, \z\ < zLAHq, of that 3D simulation, we 
have to extrapolate Vg and pg to |z| > 4i7o- To do so, we 
derive an analytic expression that fits to the Vg profile of 
the 3D simulation: 

Vg = 50cs l^tanh • (12) 

By this fitting formula, we bound the flow velocity < 
lOOcg in the high-altitude region in order to avoid unre¬ 
alistically high-speed flows. 

As for the density distribution, we adopt the following 
procedure. In \z\ < 2Ho, we use the hydrostatic density 
structure, pg = pg^mid exp(—2;^/i7g), which well explains 
the 3D MHD simulation. In \z\ > 5Ho, the 3D simula¬ 
tion shows that the mass flux PgVg = 4 x I0“®Pg_i„idCs is 
constant with 2 :. Then, we derive pg from Equation (12) 
to give the constant PgVg. We linearly connect these two 
regions between 2 ; = ±2Ho and 2 = ±3i7o- Reflecting 
that Vg approaches lOOcg in the upper regions, \z\ > 6 Hq, 
the density is also bounded by Pg > 4 x 10~^Pg^mid in 
these regions. Figure I presents the derived density and 
velocity of the gas component. 

3. RESULTS 
3.1. Time evolution 

Figure 2 demonstrates the time evolution of dust den¬ 
sity, pd, for a single grain size, a = 40pm, at r = 1 AU 
of the MMSN, or rs,mid = 7.2 x 10“® The figure clearly 
exhibits the dust grains initially located at the midplane 
are stirred upward with time. After approximately three 
rotations (red line), the dust grains are diffused by the 
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Fig. 1.— Vertical distribution of density (pg; top) and velocity 
normalized by the sound speed (ug/cg; bottom) of the background 
gas component. 

turbulence up to 2 = ±2iJo where the disk winds launch. 
Around t ^ 15 rotations (blue line), the local peaks of 
the dust-to-gas ratio gradually form at 2 ~ ±4i7o- The 
locations of the peaks approximately correspond to the 
sonic points of the background gaseous disk wind where 
the drag force switches from the small Av regime to the 
large Av regime (Equation 2; see Section 3.2 for more 
details). The location of the Pd/Pg peaks depends sen¬ 
sitively on the size of dust grains. Although in this pa¬ 
per we assume the spherical shape and the uniform solid 
density for grains, it also depends on the properties of 
dust in realistic situations. The dust-to-gas ratio further 
increases with time, and it saturates after t = 50 rota¬ 
tions (green line) to be in the quasi-steady state. The 
region inside the Pd/Pg peaks, the dust-to-gas ratio ap¬ 
proaches the volumetric average value = 0.01 because of 
the well coupling between the gas and the dust. We con¬ 
tinue the simulation further up to 200 rotations to stndy 
the quasi-steady-state behavior, which we discuss in the 
subsequent sections. 

3.2. Dependence on Grain Size 

We examine the dependence of the vertical dust distri¬ 
bution on grain size, a, or stopping time, rs,mid. Figure 3 
presents the time-averaged steady-state vertical profile of 
the dust-to-gas ratio (top panels) and the mean flow ve¬ 
locity of the dust, Vd, (Equation (11); bottom panels) for 
different stopping times. The left panels show the cases 
with Ts^mid = 1-8 X 10“®, 1.8 X 10“^, 1.8 X 10“®, and 
1.8 X 10“®, which respectively correspond to the grain 
size of a =0.1, 1, 10, and 100 pm at r =1 AU of the 
MMSN. 

In the region around the midplane, \z\ < 2Hq, the 
dust-to-gas ratio is consistent with the volumetric av¬ 
erage value, Pd/Pg = 0.01, which is independent from 
Ts.mid, because the dust is strongly coupled with the gas. 
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Fig. 2.— Time evolution of the dust-to-gas ratio, pd!Pg, for tho 
dust grains of uniform size, a = 40/im, at 1 AU. (Stopping time, 
■Ts mid = 7.2 X 10“®) Note that the value of pd/pg is multiplied by 
100. The red, blue, and green lines correspond to rotation time 
= 3, 15, and 50, respectively. After t>50 rotations, the system is 
in the quasi-steady state. 


However, the value of Pd/pg depends largely on rs,mid in 
the high-latitude regions, since the coupling is weaken 
(the local Ts increases) owing to the decrease of the gas 
density with elevating height. The dust with the small 
grain size, rs,mid < 10“®, is well coupled to the gas in 
the entire simulation domain, and is blown upward with 
the vertical outflows of the gas. On the other hand, the 
dust with the large grain size, Ts^mid > 10“®, remains in 
the disk. The velocity of the mean dust flow, Vd, in the 
wind region is faster for a smaller grain size (smaller Tg) 
owing to the well coupling to the gas. In the case with 
the largest grain size, Tg^mid = 1-8 x 10“^, Vd is directed 
to the midplane, which is however compensated for by 
the diffusion, J/pd in Equation (7). Small grains escape 
from the simulation domain, which leads to the loss of 
the total mass of the dust. However, the mass loss is 
almost negligible during the simulation time, t < 200 ro¬ 
tations, and then, it is reasonable to treat the system is 
in the quasi-steady state even for the well-coupled small 
grains. 

Among the four cases, the dust-to-gas ratio of the case 
with Tg^niid = 1.8 X 10“® shows peculiar behavior; Pd/Pg 
increases with an altitude in the high-altitude layers, 
2-5i7o ^ l-^l ^ 4.5i7o- On the other hand, the dust 
remains in the disk and pd = 0 in the high-altitude 
layers of the case with Tg^mid = 1.8x10“®. In the 
right panels of Figure 3, we focus on grain sizes with 
Tg^^^id = 4.5 X 10“®, 7.2 X 10“®, 8.1 x 10“®, which lie 
between these two cases. One can see the peaks of 
Pd/pg near z ~ ±4i7o in these cases (right top panel), 
which indicates that the dust grains float up to the lo¬ 
cation around the peaks and stagnate there. The emer¬ 
gence of the peaks can be explained for dust grains of 
Tg = 7.2 X 10“® in the following. 

The inset of the bottom right panel of Figure 3 zooms 
in the dust velocity, Vd, near the location of the peak in 
the positive z side. Vd slightly decreases from z « 3.7i7o 
to Ri 4i7o, and shows a local minimum at z Rs 4iJo, which 
almost coincides with the sonic point at z = 4.2i7o of the 
gaseous disk wind. A quasi-steady state is achieved in the 
simulation, and thus, pdVd ~ const. Therefore, the local 
minimum of Vd leads to a local maximum of pd; dust 


grains stagnate around this location. 

The non-monotonic behavior of Vd can be explained 
by the variation of the upward drag force against the 
downward gravity. Below the location of the Pd/Pg peak 
at |z| R:! 4iJo, the velocity of the gaseous disk wind is still 
subsonic, Vg < Cg. The relative velocity. At), of the dust 
to the gas is also < Cg, and the drag force is in the small 
Az; regime in Equation (2). Therefore, /d ~ 1 (Equation 
4), and Tg increases with height because of the decrease 
of pg: 




Pg{z) 


.7.2 X 10- 


\ / Pg(^)/Pg,mid \ ^ 

y ( 4 X 10“® ) 

(13) 


where we normalize with the gas density, pg = 4 x 
10“®pg_niid at the location, |z| ^ 4i7o, of the local mini¬ 
mum Vd. In |z| < 477 oj the coupling between the gas and 
the dust is lessened with an increasing height. Therefore, 
the downward gravity gradually dominates the upward 
gas drag force, which causes the deceleration of Vd in 
3.7Ho < |z| < 4Ho. 

However, Vd rapidly increases beyond |z| > 4i7o. The 
relative velocity, Av, of the dust to the gas exceeds Cg be¬ 
yond the sonic point of Vg located at |z| = A.2Hq. Then, 
the drag force switches to the large Av regime in Equa¬ 
tion (2). In the large Av :$> Cg limit, /d oc A4 (Equa¬ 
tion 4), which gives Tg oc (pgZ^th/d)”^ oc (pgAz;)“^ ^ 
{pgVg)~^ = const. (Equation 5), where we here assume 
Ad = |Dg — Dd| Vg. This shows that Tg, and therefore 
istop, does not increase (the coupling between gas and 
dust is not weaken), with an increasing altitude in the 
large Ad regime. As a result, the upward gas drag force 
increases with the increasing Vg (the 2nd term of Equa¬ 
tion 7) and dominates the downward gravity (the 1st 
term of Equation 7, which gives the rapid acceleration of 
Vd in \z\ > iHo. 

As explained above, the local peaks of the Pd/Pg ratio 
emerges when dust grains satisfy the condition that they 
are loosely coupled (ts ^ 0.1 — 0.2) to the background 
gas near the sonic point of the disk wind. By calculating 
many cases with different stopping times, the range of 
Ts.mid that gives Pd/Pg peaks is 


Ts.n.id^(4.5-8.1) X 10“®. (14) 

Smaller grains escape from the simulation domain, be¬ 
cause they are well coupled to the gas flow even near the 
sonic point. Larger grains remain in the subsonic region 
of the disk wind because they cannot be lifted up beyond 
the sonic point where the effective drag force in the large 
Ad regime could act. Therefore, dust grains outside the 
range of Equation (14) do not show local Pd/Pg peaks. 


3.3. Radial dependence 

The result based on the local vertical box (Figure 3) 
has shown that the behavior of the dust is determined 
by the nondimensional stopping time, Tg^mid, at the mid¬ 
plane. Therefore, once we specify the radial dependence, 
we can determine the dust properties in an r — z plane. 
For illustrative purpose, we adopt the MMSN (see Equa¬ 
tion 6 and below), and show the basic character of the 
vertical motions of dust grains as a function of the size a 
and the distance from the central star in Figure 4. The 
pink hatched region in the top panel indicates the radial 
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Fig. 3.— Time-averaged vertical structure of the dust distribution with various Tg j^id during t = 100 — 200 rotations. The top panels 
show the dust-to-gas ratio, and the bottom panels show the mean velocity (Equation 10). The inset in the bottom panels zooms in 
the region near = 0. The left panels present the results with Tg mid = 1-8 X 10 ® (red), 1.8 X 10 (blue), 1.8 X 10 ® (green), 
and 1.8 X 10“®(pink), which correspond to a =0.1, 1, 10, and 100 at 1 AU of the MMSN. The right panels present the cases with 
"rs.mid = 4.5 X 10“®, 7.2 X 10“®, 8.1 X 10“®, which correspond to a =25, 40, 45 /rm and show the Pd/pg peaks in \z\ ~ 4Hq. 


dependence of the size, a, of the dust grains that floats in 
the high-altitude layers with the Pd/Pg peaks. Reflecting 
the dependence of Ts^mid (Equation 6), the size of 
the floating dust grains is smaller for larger distance, 

afloat - (25 - 45) pva (y^J • (1®) 

The dust grains above the pink hatched region basically 
remain in the disk, while those below the pink hatch are 
blown upward with the vertical gas flow. 

The bottom panel of Figure 4 illustrates the pink 
hatched region in the top panel in an r — z plane. This 
panel clearly shows that dust grains are floating in a nar¬ 
row range near or slightly below the sonic point of the 
background gaseous disk wind. The size of the floating 
grains decreases for larger r. Smaller-size grains remain 
only in an outer region of the disk, while larger-size grains 
are distributed in a broader region. For instance, grains 
with a = 0.1/im remain in r > 50 AU, while grains with 
a = 100/rm remain in r > 0.5 AU. If we consider the 
time-evolution of a protoplanetary disk, the dispersal of 
the dust begins with small-size grains 
If we focus on dust grains with a certain size, say a = 
I/tm, they are blown away by the vertical gas outflow in 
the inner region, r < 10 AU, and remain in the disk in 
the outer region, r > 10 AU. Therefore, we expect that 
the evacuation of dust grains is expected to occur in an 
inside-out manner, which is qualitatively similar to the 
evolution of the gas component driven by the disk wind 
(Suzuki et al. 2010). 

4. DISCUSSION 

4.1. Effect of Reduced Mass Flux of the Disk Wind 


In this paper we use a simple model that adopts the 
time- and horizontal-plane-averaged gas density and ve¬ 
locity derived from the local 3D shearing box simulation 
solving ideal MHD equations (Suzuki & Inutsuka 2009) 
for the background flow. 

However, there is uncertainty in the mass flux of disk 
winds calculated by shearing box simulations. In partic¬ 
ular, the outflow mass flux shows an decreasing trend on 
an increasing vertical box size (Suzuki et al. 2010; Fro- 
mang et al. 2013), which implies that the mass flux of 
the disk wind we have used so far might overestimates 
an actual value. Therefore, we here check how the dis¬ 
tribution of the dust is modified when the outflow mass 
flux of the gas is smaller. In comparison to the original 
setup, which we call “Model 1” hereafter, we adopt 10 
times smaller vertical mass flux of the gas component, 
PgVg, in “Model 2”, by reducing pg by 1/10 in the wind 
region (Figure 5) but leaving Vg unchanged. 

Figure 6 presents the time-averaged steady-state ver¬ 
tical profile of the dust-to-gas ratio (top panel) and the 
mean flow velocity of the dust, Vd, (bottom panel) of 
Model 2 in comparison to those of Model 1. The gas 
density of Model 2 is 1/10 of that of Model 1 in the super¬ 
sonic region. Since Tg oc 1/pg (Equation 5), it is expected 
that in Model 2 dust grains with smaller Tg^mid at the 
midplane will give similar results to Model 1 in the su¬ 
personic region. Therefore, we adopt Tg^mid = 7-2 x 10“^ 
in Model 2, which is 1/10 of Tg^i„id(= 7.2 x 10“®) of Model 
1, whereas Tg^mid = 7.2 x 10“® and 7.2 x 10“^ correspond 
to the grain size of a = 40/tm and a = 4/rm at r = 1 AU 
of the MMSN, respectively. 

Figure 6 shows that the velocity profiles of these two 
cases are almost the same, and so are the locations of 
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Fig. 4.— Properties of the dust grain dynamics in the MMSN 
with the disk wind. Top: Behavior of the dust with different grain 
size, a, on radius from the central star. The pink hatched region 
corresponds to the dust particles that float in the high-altitude 
layer and show a sharp concentration of pd/pg around l^l ~ 4//o- 
The dust particles below the pink hatched region are blown away 
with the gas outflows from the simulation domain. Those between 
the pink hatch and the green line are distributed in the disk without 
pd/pg peaks but are stirred up to \z\ > IHo. Those above the green 
solid line are settled within IHq around the midplane. Bottom: 
Floating dust grains in an r — 2 diagram. Color contour indicates 
the grain size in the pink hatched region in the top panel. The 
black solid line indicates the sonic point of the gaseous disk wind, 
and the dotted line denotes one scale height, 2 : = IHq. 



Fig. 5.— Comparison of the gas densities of Models 1 and 2. 

the peaks. This indicates that the dust floating by the 
stagnation of upflowing grains is a robust feature in a 


qualitative sense even though the mass flux of the gaseous 
disk wind might be reduced, whereas the size of floating 
dust grains is smaller for smaller mass flux of the disk 
wind. Model 2 shows larger pdjpg in the subsonic region, 
1^1 < 4:Ho, because the gas density of Model 2 there is 
not as low as 1/10 of that of Model 1. Thus, Ts in the 
subsonic region of Model 2 is smaller than Ts of Model 1, 
and a larger amount of dust grains are lifted up by more 
efficient coupling to the gas. 



Fig. 6.— The same as the right panels of Figure 3 but for the 
comparison between Models 1 and 2. 


4.2. Limitations 

In the previous subsection, we have discussed the un¬ 
certainty of the wind mass flux in the local shearing box 
approximation. In addition to this, the configuration of 
large-scale magnetic field is another effect which cannot 
be taken into account in the local treatment. Gas flows 
along large-scale magnetic field lines, and the mass con¬ 
servation should follow pgVg/Bp = const, in realistic sit¬ 
uations, where Bp is poloidal magnetic field. Thus, the 
mass flux, PgVg, would probably decreases as the mag¬ 
netic flux tube expands with altitude, and hence Bp de¬ 
creases. As a result, the location and detailed shape of 
PdlPg peaks might be modified in such situations. How¬ 
ever, the physical mechanism that forms a Pd/Pg peak 
near the sonic point of vertical npflows (Section 3) is ro¬ 
bust in a qualitative sense. 

In this paper, we have assumed the steady-state con¬ 
dition for the background gaseous flows. In reality, how¬ 
ever, the disk wind is intermittent with quasi-periodicity 
of 5-10 rotations (Suzuki & Inutsuka 2009, 2014). We 
need to take into account effects of time-dependent ver¬ 
tical upflows on dust grains in more detailed studies 
(Turner et al. 2010). 

The strong coupling approximation adopted in this pa¬ 
per (Equation 9) is basically acceptable for floating dust 
grains because Ts{~ 0.2) is kept less than unity near and 
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above the location where they are floating (Section 3). 
However, if time-dependent vertical upflows were taken 
into account, grains with Tg > 1 might appear locally in 
a transient manner. For those grains, it is necessary to 
solve the time-dependent equation of motion (Equation 

7 ). 

We treat dust grains as test particles, because we start 
the simulations from the total dust-to-gas ratio = 0.01, 
and the condition, Pd/Pg < 1, is satisfied in the through¬ 
out the simulations. However, when we consider the de¬ 
pletion of the gas component at later evolutionary stages, 
Pd/pg might exceeds unity especially near the peak loca¬ 
tion. In order to treat the long-time evolution of proto¬ 
planetary disks, the backreaction from the dust on the 
gas (e.g., Johansen & Klahr 2005) 

We use the result of the ideal MHD simulation, namely 
we do not consider the effect of an MRI-inactive re¬ 
gion (“dead zone”; e.g., Gammie 1996). If the effects 
of non-ideal MHD are taken into account, the turbu¬ 
lence strength, a, would be smaller than the adopted 
value, a = 0.01, and thus it takes time for grains at 
the midplane at the beginning of the calculations to dif¬ 
fuse upward (see Section 2.1). After the quasi-steady 
state is achieved, however, a possible existence of dead 
zone would not change our result on floating dust grains 
at least in a qualitative sense, because the disk wind is 
driven from the MRI-active region near the surface where 
the sufficient ionization is attained. For instance, Suzuki 
et al. (2010) demonstrated that the mass flux of the disk 
wind is weaken by half for a typical case with a dead 
zone. Recent simulations that include various non-ideal 
MHD effects also show that magnetically driven winds 
are launched in spite of the strong magnetic diffusion 
(e.g., Bai & Stone 2013b; Gressel et al. 2015). Although 
the size of floating grains might be smaller than that 
predicted from our results based on the ideal MHD sim¬ 
ulations, floating dust grains would be a robust feature 
even if non-ideal MHD effects were considered. 

When we consider dust grains with a = 0.1 pm at r = 1 
AU, the dispersal timescale is « 7000 yrs, which is much 
shorter than the typical lifetime of protoplanetary disks, 
« 106-7 yj.g (Haisch et al. 2001). The main reason of 
this discrepancy is that the effect of the radial accre¬ 
tion cannot be taken into account in the local box. In 
realistic situations, dust grains would be supplied from 
the outer region and the actual dispersal time would be 
much longer. Suzuki et al. (2010) have done such similar 
calculations for the gas component and showed that the 
dispersal of the disk gradually takes place with timescale 
of ~ 10® yrs even though the dispersal time at r = 1 
AU is Rs 4000 yrs in the local treatment. The similar 
argument is applicable to the dust component, and we 
expect that the dispersal timescale of small dust grains 
will turn out to be a reasonable value in a global model. 

Since we have restricted our calculations to those in 
the local vertical box in this paper, local concentrations 
of floating dust grains are formed by the combination 
of the vertical component of the gravity and the drag 
force. In realistic situations, however, the direction of 
the flow, the gravity, and the drag force are not gen¬ 
erally coaligned each other. In such circumstances, up- 
flowing grains could play a role in large-scale circulation 
of the dust in protoplanetary disks. The disk wind of 
the gas component is probably accelerated by magneto¬ 


centrifugal force to the radially outward direction (e.g., 
Blandford & Payne 1982). Dust grains that are well cou¬ 
pled to the gas are uplifted and transported to larger r 
along large-scale magnetic field by such gaseous upward 
and outward flows. Eventually, the ambient gas den¬ 
sity decreases as the magnetic flux tube expands with 
altitude, which makes the coupling to the gas weaken 
to give Ts > 1. Those grains will settle down to the 
midplane by the vertical component of the central-star 
gravity without sufficient gas drag force in such a distant 
region. Outflows are expected to blown out intermit¬ 
tently (Suzuki & Inutsuka 2009), and hence, dust grains 
are lifted up and settled down in a time-dependent man¬ 
ner in realistic situations, rather than the formation of 
time-steady local peaks discussed in the previous sec¬ 
tions. The dust grains settled down at larger r further 
accrete inward, being coupled to the accreting gas. Dust 
grains are also supposed to grow through this large-scale 
circulation (Shu et al. 1996). 

4.3. Observational Implications 

Finally we discuss observational implications of float¬ 
ing dust in protoplanetary disks. McJunkin et al. 
(2014) estimated the column densities of neutral hy¬ 
drogen around many T-Tauri Stars and Herbig Ae/Be 
stars by the circumstellar adsorptions of the Lyman-a 
emissions from the central stars. The column densities 
derived from Lyman-a absorption are smaller than the 
values derived from the dust extinction measurements. 
A possibility of a high dust-to-gas ratio in the upper- 
atmosphere of the protoplanetary disks was proposed as 
a possible solution to the discrepancies. Ellerbroek et al. 
(2014) observed outflows/winds from a Herbig Ae star. 
The object shows variability where optical fading events 
are correlated with near-infrared brightening, which im¬ 
plies that a sizable amount of the dust grains are con¬ 
tained in the variable outflows/winds. The mechanism 
studied in this paper may provide an explanation for ex¬ 
istence of those dusty materials in upper atmospheres of 
protoplanetary disks. 

Also, the scattered light image in circumstellar envi¬ 
ronment of RY Tau by near-infrared coronagraphic imag¬ 
ing polarimetry indicated that the observed polarized 
intensity distribution shows a butterfly-like distribution 
as shown by the observation at millimeter wavelengths 
(Takami et al. 2013). They performed comparison be¬ 
tween the observed polarized intensity distribution and 
models consisting of gaseous disk and dust. The result 
suggests the scattered light in the near-infrared is associ¬ 
ated with an optically thin and geometrically thick layer 
at an upper layer. This phenomenon could be explained 
by the dust grains lifted up by the disk wind. In our 
future studies, we plan to create pseudo observational 
images (e.g., Takami et al. 2014) based on a theoretical 
model for floating dust grains. 

5. SUMMARY 

We have studied the dynamics of dust grains in the 
protoplanetary disk winds driven by the MRI-triggered 
turbulence. Our calculation in the ID vertical box 
has shown that the intermediate-size dust grains with 
Ts.mid 7.2 X 10“® float at about 4 scale heights from the 
midplane. This location corresponds to the sonic point of 
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the background gaseous disk wind, where the drag force 
from the gas to the dust changes from the small relative 
velocity regime to the large relative velocity regime. The 
size of the floating dust grains depends on the distance 
from the central star; if we take the MMSN for exam¬ 
ple, we have size of floating dust, afloat oc smaller 

grains float at larger r. The dust grains with the smaller 
size than that of the floating ones are well coupled to the 
bulk gas flow and continue to stream out of the simula¬ 
tion domain, while the larger dust grains remain in the 
disk; this mechanism selectively disperses small grains in 
an inside-out manner and leaves large grains in proto¬ 
planetary disks. 

One of the interesting findings in this paper is the char¬ 
acteristic size of the dust grains floating and stagnating 
at an upper atmosphere of the protoplanetary disk. The 
size range of those grains is a function of the distance 
from the central star, i.e., pg{r),Cs{r) and appears to be 


relatively narrow as shown in the top panel of Figure 4. 
Future observations of the emission and/or absorption 
by those floating dust grains may provide information 
about the typical size of those grains; i.e., the detection 
of a dust floating layer may immediately enable us to de¬ 
termine the size of those particles. Combining with the 
possible observations of the disk temperature , i.e., Cs(r) 
and density, we will be able to constrain the mass flux of 
bulk flow, PgWg, which is difficult to observe directly. This 
may provide important observational diagnostics on the 
dispersal timescale of individual protoplanetary disks. 
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